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ABSTRACT 

We report on the results of R band observations of the error box of the 7-ray 
burst of August 28, 1997, made between 4 hours and 8 days after this burst 
occurred. No counterpart was found varying by more than 0.2 magnitudes down 
to R — 23.8. We discuss the consequences of this non-detection for relativistic 
blast wave models of 7-ray bursts, and the possible effect of redshift on the 
relation between optical absorption and the low-energy cut off in the X-ray 
afterglow spectrum. 

Subject headings: Gamma-rays:bursts — gamma-rays:observations— 
galaxies: general — radiation mechanisms:non-thermal — dust, extinction 
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1. Introduction 

Since the first discovery of a 7-ray burst (GRB) in 1967 (Klebesabel et al., 1973) 
these short outbursts of highly energetic photons have formed one of astronomy's most 
elusive problems. Following the discovery by Meegan et al. (1992) of their isotropic sky 
distribution and inhomogeneous spatial distribution (which excluded that GRBs originate 
from a galactic-disk source population) the discussion on the nature of GRB sources 
focussed on their distances: either of order 10 5 pc ('galactic halo' model), or several Gpc 
('cosmological' model). The association of the optical counterpart of GRB970228 (Groot et 
al., 1997; Van Paradijs et al., 1997) with what is most likely a galaxy (Groot et al, 1997b; 
Metzger et al., 1997a; Sahu et al., 1997) and especially the determination of a redshift for 
GRB970508 (Metzger et al., 1997b) have shown that GRBs are located at cosmological 
distances, and are thereby the most luminous photon sources known in the Universe. The 
question of what causes GRBs has now become the centerpiece of the discussion, and the 
detection of more optical counterparts is a key element in determining this cause. 

In this Letter we report on our search for a transient optical counterpart for 
GRB970828, based on observations made with the 4.2m William Herschel Telescope (WHT) 
on La Palma, and the 3.5m WIYN Telescope on Kitt Peak. None was detected, down to a 
magnitude level R = 23.8. 

GRB 970828 was discovered with the All-Sky Monitor (ASM) on the Rossi X-ray 
Timing Explorer (RXTE) on August 28, 1997, UT 17 h 44 m 36 s from an elliptical region 
centered at RA=18 h 08 m 29 s , Dec=+59°18^0 (J2000), with a major axis of 5f0, and a minor 
axis of 2:0 (Remillard et al. 1997; Smith et al. 1997). Within 3.6 hours the RXTE/PCA 
scanned the region of the sky around the error box of the ASM burst, and detected a weak 
X-ray source, located in the ASM error box with a 2-10 keV flux of 0.5 mCrab (Marshall 
et al. 1997). The burst was also detected with the Burst And Transient Source Experiment 
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BATSE and the GRB experiment on Ulysses. Its fruence and peak flux were 7 x 1CT 5 
erg cm' 2 , and 3 x 10~ 6 erg cm -2 s" 1 , respectively. From the difference between burst 
arrival times, its position was constrained to lie within a 1.62 arcminute wide annulus, 
that intersected the RXTE error box (Hurley et al. 1997). In an ASCA observation made 
between Aug. 29.91 and 30.85 UT, a weak X-ray source was detected at an average flux level 
of 4 x 1(T 13 erg cm' 2 s' 1 (2-10 keV). The ASCA error box is centered on RA=18 h 08 m 32!3, 
Dec =+59°18'54"(J2000) and has a 0^5 radius (Murakami et al. 1997). 



2. Observations and Data Analysis 

We observed the GRB error box with the Prime Focus Camera of the WHT, on 9 
nights between August 28, UT 21 h 47 m , and September 5, UT 22 h 07 m (see Table 1). The 
first observation was made just over 4 hours after the 7-ray burst. All observations were 
made with a Cousins R band filter (Bessell 1979). During the first two nights and the last 
three nights, we used a LORAL 2048x2048 CCD chip, with 15/x pixels, giving a field of 
view of 8'45x8f45. During the intervening nights we used an EEV CCD chip (2048x4096), 
windowed at 2048x2400, with 13.5/x pixels giving a 8'lx9f5 field of view. On August 30 two 
R band images were made with the WIYN Telescope. The camera contained a 2048 x 2048 
CCD, giving a field of view of 6^8 x 6f8. 

We obtained a photometric calibration of the CCD images from observations of Landolt 
Selected Area 113, stars 281, 158, 183 and 167 (Landolt 1992), on Aug 31, h 14 m UT with 
the WHT. 

A region of 2'x2' centered on the ASCA position in the bias-subtracted and flatfielded 
images was analyzed using DoPhot (Schechter et al., 1993), in which astrometric and 
photometric information of all objects are determined from bivariate Gaussian function 
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fits to the brightness distribution in their image; the parameters of these fits also tell us 
whether an object is stellar (i.e., unresolved) or a galaxy. In this region (see Fig. 1) we find 
a total of 63 objects, 36 of which are stellar, and 27 galaxies, down to i?=23.8. 

We have searched for variable objects by comparing the magnitudes of each star as 
determined for each of the images. Comparison of images taken on different nights showed 
no variation on time scales between a day and a week in excess of 0.2 mag for R <23.8 (for 
the last three nights the limit on variability is 0.3 mag for R <23.8). Comparison of three 
images taken on the night of August 29 to 30 showed no variations on time scales of several 
hours in excess of 0.2 mag for R < 22.5. 

3. Discussion 

3.1. Comparison with optical afterglows of GRB970228 and GRB970508 

The large variation in optical response of GRBs (relative to their strength in 7 rays) 
was already clear from a comparison between GRB970228 and GRB970111. Within a day 
after GRB970228 ocurred it showed an optical afterglow at R = 20.8 (Van Paradijs et 
al. 1997; Galama et al. 1997a; Pedichini et al. 1997; Guarneri et al. 1997). GRB970111 
was not detected in optical observations made 19 hours after it occurred (R > 20.8, and 
R > 22.6, for variations in excess of 0.2 and 0.5 magnitudes, respectively, Castro-Tirado et 
al. 1997), in spite of the fact that its 7-ray fluence (Galama et al., 1997b) was five times 
larger than that of GRB970228 (Costa et al., 1997). Since only one deep image was made 
in the week following GRB970111, its non-detection may have been the result of, e.g., a 
very rapid decay of any optical afterglow, or a very slow rise thereof (like for GRB970508, 
see Bond 1997; Djorgovski et al. 1997; Sahu et al. 1997b; Galama et al. 1997c). 

The non-detection (R > 23.8 for variations in excess of 0.2 magnitudes) of GRB970828 
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during our optical observations, which covered the time interval between 4 hours and 8 
days after the burst at intervals of a day, show the very large range in optical responses of 
GRBs in an even more striking fashion. We have used the fluence, -Eqrb (in ergs cm -2 ), 
as a measure of the GRB strength, and compared the ratio of the optical peak flux to the 
GRB fluence of GRB970828 with that of GRB970508. The latter had a peak magnitude 
R = 19.8 (Mignoli et al., 1997), therefore the difference in optical peak luminosities between 
GRB970508 and GRB970228 is more than 4 magnitudes. The ratio of their fluences, 
Eqrb (970828) /E GRB (970508) = 24 (Kouveliotou et al. 1997a,b). Thus, we find that the 
optical peak response of GRB970828, with respect to its 7-ray fluence, is a factor ~ 10 3 
smaller than that of GRB970508. (Compared to GRB970228 the difference is a factor 
> 10 2 .) 

We have made a similar comparison with published X-ray afterglow fluxes (F x ) for the 
two GRBs with optical afterglow. Most of these refer to the energy range 2-10 keV. Only 
the ROSAT fluxes had to be transformed to this range; in doing this we assumed a power 
law X-ray spectrum with photon index in the range -1.4 to -2.0 (Costa et al. 1997; Yoshida 
et al. 1997). This range leads to an uncertainty in the transformed ROSAT flux of less than 
a factor 2. The results, in the form of the ratio Tlx = Fx/E GRB , are summarized in Fig. 
2, which shows the variation of this quantity as a function of the time interval since the 
burst, for four bursts with published X-ray afterglow information. This figure shows that 
the differences in Tlx between these bursts are moderate (less than an order of magnitude). 
It is noteworthy that the two bursts with optical counterparts also have the highest values 
of Tlx (for a given value of At). 

We finally compared the peak flux in the i?-band afterglows with the brightness of 
the X-ray afterglow. In view of their rather similar decay rates we used for the latter the 
2-10 keV flux as measured 1 day after the GRB occurred, i*x(l day). The corresponding 
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ratio F pcak ( J Rband)/F x (l day) for GRB970828 differs by a factor > 150 from that for 
GRB970508, and a factor > 10 from that for GRB970228. 

3.2. Comparison with Relativistic Blast Wave Models 

A relatively succesful way of explaining the existence of GRB afterglows (at all 
wavelengths) has been the so-called blastwave or fireball models (e.g. Meszaros 1995). 
These models involve the generation of a massive amount of energy in a very small, 
compact region, by an unexplained mechanism. The result of this dumping of energy 
is a relativistically expanding fireball (blastwave), that collides with the interstellar or 
circumstellar medium and generates shocks that emit the synchrotron radiation that is 
observed as the afterglow. 

Figure 3 shows the available data for GRB970828 in 7-rays, X-rays, B, and R, plus 
simple blastwave model fits, which are normalized to agree with the X-ray data. If we 
compare this with the data available for GRB970228 (Wijers , Rees and Meszaros, 1997), it 
is striking that the decay part of the X-ray curves are virtually the same for these two bursts 
(i.e. in slope and offset). But whereas the first stages of the optical decay for GRB 970228 
are in good agreement with the afterglow prediction (Wijers, Rees, and Meszaros 1997), 
the earliest upper limit to the optical brightness of GRB 970828 is 300 times lower than the 
predicted value. 

The simplest spherically symmetric blastwave models for GRB afterglows require that 
the slope of the spectrum follows from the slope of the temporal decay, once the decay curve 
is measured in one wavelength band and is found to be a pure power law. From that, the 
offset in brightness at any other waveband is fixed and the predicted flux at that waveband 
is hard to change. Meszaros, Rees, & Wijers (1997) showed that if the blastwave is beamed 
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one can get different relations between spectral and temporal slopes, giving possibly much 
smaller offsets between the optical and X-ray light curves of the afterglow. As an example, 
let the energy per unit solid angle, E, vary with angle from the jet axis, 9, as E oc 6~ A and 
the Lorentz factor T oc 9" 1 . Then a temporal decay rate F oc t~ 13 as seen here would occur 
for a spectrum F v oc u 0A , i.e. it would rise from optical to X rays and the predicted R band 
curve would be a factor 24 below the X-ray curve. At the time of our first limit this model 
would give R = 28.4, quite consistent with the data. 



3.3. Absorption in Redshifted Material 

Another explanation, pointed out to us by dr. B. Paczyhski, for the non-detection of 
optical afterglow could be photoelectric absorption, also visible as a low-energy cut-off in 
the X-ray spectrum. If we assume a modest hydrogen column density of Ah ~ 10 21 atoms 
cm" 2 and make the assumption that the absorbing material is at redshift z—0, this would 
imply 0.34 magnitudes of extinction in the R band (Gorenstein 1975; Cardelli et al. 1989). 

In case the absorption takes place at some redshift z the effect is a bit more complicated. 
The cross section for photo-electric absorption in the (0.2-5) keV range depends on energy 
roughly as E~ 2S (Morrison and McCammon 1983). Then the factor by which the apparent 
Ah, inferred from the low-energy cut-off in the X-ray spectrum, has to be increased is 
approximately (1 + z) 2S . If we assume, for example, that the GRB occured at a redshift 
of z=l, the factor by which the apparent value of A H has to be increased would be ~6. 
Moreover, the photons in the R band we observe would be at wavelengths near 3200 A at 
the source, at which wavelength the interstellar absorption is approximately a factor 2.5 
larger than in the R band (Cardelli et al., 1989). These combined effects would lead, for a 
GRB at z=l and an apparent, moderate, Ah=10 21 atoms cm -2 to an R band extinction of 
~5 mags. 
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If absorption is the correct explanation, a substantial fraction of GRB sources (those 
with a very small optical response) would be located close to where large column densities 
are available, i.e., in disks of galaxies. This would link GRBs to a population of massive 
stars. This is expected for the failed-supernova model and for the hypernova model, 
proposed by Woosley (1993) and Paczyhski (1997), respectively. In view of the large kick 
velocities imparted on neutron stars at birth (Lyne & Lorimer 1994; Hansen & Phinney 
1997; Van den Heuvel & Van Paradijs 1997) it remains to be seen whether a merging 
neutron star binary model would be consistent with this consequence. 
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Figure Captions 

Fig. 1: 2'x2'i?-band image of the sky region centered on the 0.5' radius ASCA error 
box of GRB970828, taken at the WHT on Sept 2. 

Fig. 2: Variation of the ratio Rx of the afterglow (2-10 keV) X-ray flux (see text) to 
the fluence in the 7-ray burst as function of time in days. GRB fluences were obtained from 
the following sources. GRB970228: Costa et al. (1997); GRB970508: Kouveliotou et al. 
(1997a); GRB970828: Kouveliotou et al. (1997b). The X-ray fluxes were obtained from the 
following sources: GRB970228: Costa et al. (1997), Yoshida et al. (1997), Frontera et al. 
(1997); GRB970508: Piro et al. (1997); GRB970828: Marshall et al. (1997), Murakami et 
al. (1997b). 

Fig. 3: Variation of the observed fluxes in 7 rays, X rays, and in the B and R bands of 
GRB970828, together with simple blast wave model fits as described in Wijers, Rees and 
Meszaros (1997). 



Table 1: Log of observations GRB970828 
Date Telescope UT Start Exp. time (s) Seeing 



Aug 28 


WHT 


21:47 


900 


0"86 


Aug 29 


WHT 


21:15 


900 


CC/74 


Aug 30 


WIYN 


05:08 


600 


0"8 


Aug 30 


WIYN 


07:38 


900 


1"2 


Aug 30 


WHT 


23:22 


900 


0'.'90 


Aug 31 


WHT 


20:54 


900 


0"71 


Sep 1 


WHT 


21:16 


600 


0"80 


Sep 2 


WHT 


20:53 


600 


0"76 


Sep 3 


WHT 


22:44 


600 


0"88 


Sep 4 


WHT 


21:53 


600 


0"79 


Sep 5 


WHT 


22:07 


600 


0"86 
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